< 



New low-frequency nonlinear electromagnetic wave in a 

magnetized plasma 

M. Marklund, P.K. Shukla, L. Stenflo, G. Brodin, and M. Servin 

Department of Physics, Umed University, SE-901 87 Umed, Sweden 

O ■ ABSTRACT 

o : 

(N ■ 

A new nonlinear electromagnetic mode in a magnetized plasma is predicted. 
Its existence depends on the interaction of an intense circularly polarized electro- 
magnetic wave with a plasma, where quantum electrodynamical photon-photon 
scattering is taken into account. This scattering gives rise to a new coupling 
^ ' between the matter and the radiation. Specifically, we consider an electron- 

>■ ' positron plasma, and show that the propagation of the new mode is admitted. It 

could be of significance in pulsar magnetospheres, and result in energy transport 
^ 1 between the pulsar poles. 

Subject headings: Plasmas — pulsars: general — stars: neutron — waves 

Astrophysical environments can be most violent and energetic. Physics considered 'ex- 
otic' in Earth based laboratory applications can be common throughout our Universe, and 

+3 • sometimes even vital for the existence of certain observed phenomena. Pulsars, surrounded 
by strong magnetic fields, are most prolific sources of exotic physics. Quantum electrodynam- 
ics (QED) is an indispensable explanatory model for much of the observed pulsar phenomena. 

^ ■ Scattering of photons off photons is predicted by QED, and it can be a prominent component 
of pulsar physics, since pulsars offer the necessary energy scales for such scattering to occur. 
Related to the scattering of photons is the concept of photon splitting in strong magnetic 
fields (Adler 1971). It has been suggested that such effects could be important in explaining 
the radio silence of magnetars (Kouveliotou 1998; Baring & Harding 2001). In the present 
Letter we will point out the existence of a new low-frequency nonlinear electromagnetic wave 
that may exist in a magnetized plasma, due to the interaction of photons with the quantum 
vacuum. A discussion of the properties of this new electromagnetic wave using parameters 
relevant to strongly magnetized pulsars will be given. 

The weak field theory of photon-photon scattering can be formulated in terms of the 
effective Lagrangian density 
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where Jzf = — \^oF ab F ab = ^e (E 2 — c 2 B 2 ) is the classical free field Lagrangian, and 



HE — ft^c 



4(\F ab F ab ) 2 + 7(lF ah F ab )' 



(2) 



is the Heisenberg-Euler correction (Heisenberg & Euler 1936; Schwinger 1951), where F ab = 
^ abcd F cd , and \F ab F ab = -cE • B. Here k = 2a 2 h 3 /45m*c 5 « 1.63 x lCr 30 ms 2 /kg, a is the 
fine-structure constant, h is the Planck constant, m e is the electron mass, and c is the speed 
of light in vacuum. With F ab = d a A b — d b A a , A b being the four-potential, we obtain, from 
the Euler-Lagrange equations, the field equations d b [d^f /dF ab ] = 0, i.e. (see, e.g. Shukla et 
al. 2004) 
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where j a is the four current. 



For a circularly polarized wave E = E^x ± iy) exp(ik • x — iut) propagating along a 
constant magnetic field B = B z, the invariants satisfy 
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where k is the wave number and uj the frequency of the circularly polarized electromagnetic 
wave. Thus, Eq. (3) can be written as 



UA a = -4e K 
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in the Lorentz gauge, and □ = d a d a . For circularly polarized electromagnetic waves propa- 
gating in a magnetized cold multicomponent plasma, the four current can be 'absorbed' in 
the wave operator on the left-hand side by the replacement 



n^-D(u,k), (6) 

where D is the plasma dispersion function, given by (see, e.g. Stenflo (1976) and Stenflo & 
Tsintsadze (1979)) 

D(lu, k) = k 2 c 2 - uo 2 + Y W(J ** , (7) 
where the sum is over the plasma particle species j, 



- 3- 



cv cj = -j— and u pj = —-L , (8) 
is the gyrofrequency and plasma frequency, respectively, and 

7; = (1 + (9) 
is the the gamma factor of species j, with Vj satisfying 

^ff^y \±a (10) 

' \cm 0j ) [w(l + ijyv ±u><#' 1 ' 

Here noj denotes particle density in the laboratory frame and rriQj particle rest mass. It should 
be emphasized that for the case of circularly polarized waves propagating along an external 
magnetic field, the effects due to the electron and ion currents can be calculated without 
linearizing the plasma governing equations (see e.g. Derby 1978, Goldstein 1978, Stenflo & 
Tsintsadze 1979 and Stenflo & Shukla 2001). Thus, in Eq. (7), the relativistic nonlinearity, 
the full Lorentz force as well as the other plasma nonlinearities are fully accounted for. The 
only limiting assumption of the electromagnetic field amplitude is due to the validity of 
the Euler-Heisenberg Lagrangian. We note that the application of Eq. (2) requires field 
strengths smaller than the Schwinger critical field E s = m 2 e c'/eh ~ 10 18 V/m (see e.g. 
Schwinger 1951). The dispersion relation, obtained from Eq. (5), then reads 
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We note that as the plasma density goes to zero, the effect due to photon-photon scattering, 
as given by the right-hand side of Eq. (11), vanishes, since then J 1 — k 2 c 2 = 0. 

It may be instructive to first consider the small amplitude limit of Eq. (11). In a 
magnetized pair-plasma this reads 
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which shows the existence of fast and slow EM waves. In the low-frequency long wavelength 
limit, and for oj p <C ou c (as we typically have in pulsars) we then have u/k ~ c. The phase 
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speed of the slow wave will thus be close to the speed of light. However, for large amplitude 
waves, the situation is drastically different. As we will see below, we indeed have then 
quite new kinds of dispersive electromagnetic waves whose frequencies depend on the wave 
amplitudes, and for which the phase velocities fulfill u/k <C c. 

Next, we focus on mode propagation in an ultra-relativistic electron-positron plasma 
(7 e 3> 1), where the two species have the same number density n . Then Eq. (11) gives 
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Following Stenflo & Tsintsadze (1979), we have defined uj e = eE /cm 0e . 



Looking for low-frequency modes, we now use the approximation uj <C kc, at which Eq. 
(13) gives 
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It is sometimes advantageous to use the relation w £ = u e (E / E s ), where u e = m e c 2 /h is the 
Compton frequency, to write Eq. (14) as 
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Using the dispersion relation (14) the group velocity v g = dui/dk is 
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where v p = uj/k is the phase velocity. 



(16) 



Pulsar magnetospheres exhibit extreme field strengths in a highly energetic pair plasma. 
Ordinary neutron stars have surface magnetic field strengths of the order of 10 6 — 10 9 T, while 
magnetars can reach 10 10 — 10 11 T (Kouveliotou 1998), coming close to, or even surpassing, 
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energy densities e Eg corresponding to the Schwinger limit. Such strong fields will make the 
vacuum fully nonlinear, due to the excitation of virtual pairs. Photon splitting can therefore 
play a significant role in these extreme systems (Harding 1991; Baring & Harding 2001). 
The emission of short wavelength photons due to the acceleration of plasma particles close 
to the polar caps results in production of electrons and positrons as the photons propagate 
through the pulsar intense magnetic field (Beskin 1993). Given the Goldreich- Julian density 
ncj = 7 x 10 15 (0.1 s/P)(S/10 8 T) m~ 3 , where P is the pulsar period and B the pulsar 
magnetic field, the pair plasma density is expected to satisfy n = Muqj, where M is the 
multiplicity (Beskin 1993; Luo et at. 2002). Moderate estimates give M = 10 (Luo et at. 
2002). Thus, the density in a pulsar pair plasma can be of the order 10 18 m~ 3 . The plasma 
experiences a relativistic factor ~ 10 2 — 10 3 (Asseo 2003). On the other hand, the primary 
beam will have n ~ n GJ and 7 ~ 10 6 — 10 7 (Asseo 2003). 

Furthermore, we note that the pulsar magnetosphere is very far from thermodynamical 
equilibrium, and consequently there exist numerous possible excitation mechanisms. In 
addition to Cherenkov excitation and cyclotron excitation that may excite Alfven wave 
like modes (e.g. Lyutikov 1999), free energy sources in the form of temperature anisotropics 
and/or electron beams can drive Weibel instabilities and beam type instabilities, respectively. 
Several of these initial excitation mechanisms can lead to the existence of arbitrarily large 
amplitude waves. 

The two QED terms in the squared bracket of Eq.(15) correspond to the wave field 
and the static background field respectively. The former dominates if the oscillating (wave) 
part of the magnetic field is larger than the static part, and the latter dominates if the 
opposite is true. Due to the dipole nature of the static field, the former case is more likely to 
hold at a comparatively large distance from the pulsar surface, whereas close to the surface 
the static part is likely to dominate. From now on we concentrate on comparatively weak 
background fields strengths (thus excluding regions in the immediate vicinity of the pulsar 
surface). Accordingly we have cBq <C E$, and we therefore drop the term proportional to Bq 
in Eq. (15). Next, using the normalized quantities Vt = cucu e /ujp e , K = (4a/4:5n)~ 1 /' 2 kcuj e /u}p e 
and £ = (4o:/4:5ii)Eq/E s , the dispersion relation (15) reads 

Q2 = £2r2 t Mi- (17) 

This dispersion relation describe three different modes, two with + polarization and one with 
— polarization. The normalized frequency as a function of K and £ is shown in Fig. 1. We 
note that for K <C 1, the dispersion relation (17) agrees with that of Stenflo &Tsintsadze 
(1979), whereas in the opposite limit K ^ 1, the QED term in (17) is dominating. For the 
given density, the latter regime applies, except for extremely long wavelengths (> 10 8 m), 
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and thus we note that QED effects are highly relevant for the propagation of these modes in 
the pulsar environment. For small K there is only one mode, but as seen from the second 
and third panels of Fig. 1, two new modes appear for K > 2.6. Thus for large K, applicable 
in the pulsar environment, there are three low-frequency modes (uj <C kc) that depend on 
nonlinear QED effects for their existence. Using cBq <C Eg, the expression (16) for the group 
velocity becomes 

dtt _ Q±2Q 2 /SK 2 D, 

dK ~ &±?>tt 2 /2£K 2 K' ^ ' 

and thus we see that the propagation speed depends nonlinearly on the plasma parameters. 
Furthermore, from (18), and recalling the normalizations used, we see that the group- velocity 
is typically well below the speed of light. We emphasize that the low group velocity by itself 
makes our mode more long-living, since the energy supplied to the mode in order to balance 
the convective energy loss is proportional to the group velocity. We thus suggest that the 
three new modes presented above can contribute to an understanding of the very complicated 
energy transport phenomena taking place in the accretion discs of pulsars. 

In summary, we have reported the existence of a new low-frequency nonlinear elec- 
tromagnetic wave in a magnetized plasma. The dispersion relation of the wave has been 
presented, and analysed using relevant astrophysical parameters. Applications to pulsar 
magnetoplasmas have been pointed out. 
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Fig. I. — Dispersion surfaces Vt = Vt(K, £) as given by Eq. (17). The first panel corresponds 
to the — sign in Eq. (17), and exists for all K and £. The second panel shows the fast + 
polarized mode, which exists for K > 2.6. The third panel depicts the slow + polarized 
mode, also for K > 2.6. 



